Spectral functions of 2D systems with coupling of electrons 
to collective or localized spin degrees of freedom 
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The spectral properties of itinerant 2D systems with (nearly) ferromagnetic ground state are 
studied within the spin-fermion and the classical s-d exchange models. While the former model 
describes the effect of collective magnetic excitations on the electronic properties, the latter one 
considers the effect of local moments. We use the equation of motion approach combined with 
the 1/M and 1/z expansions (M is the number of spin components, z is the coordination num- 
ber) to investigate spectral functions. In both the models the spectrum splitting occurs at low 
temperatures T in the renormalized classical regime due to strong magnetic fluctuations. For the 
interaction J between electronic and magnetic degrees of freedom comparable to the bandwidth 
(the intermediate-coupling regime) the 1/z expansion predicts full splitting of electronic spectrum 
below some temperature T s . At the same time, the spectrum remains only partly split in the 1/M 
expansion with the quasiparticle structure which becomes fully coherent at low T. 



I. INTRODUCTION 

The narrow-band d-systems pose the problem of mag- 
netism of electronic systems, where the electron-electron 
interactions are comparable to the bandwidth. Two dif- 
ferent points of view on magnetic excitations in such sys- 
tems can be considered. On one hand, magnetic exci- 
tations can be treated as collective excitations of itin- 
erant electrons [lj]. On the other hand, in many d- and 
especially f-electron systems one can consider localized 
(magnetic) and itinerant excitations as two independent 
degrees of freedom, interacting with each other. 

The latter viewpoint is convenient, e.g., for tran- 
sition metals containing both conducting and local- 
ized electrons. For the description of such systems 
the s-d exchange model was introduced about 50 years 
ago by Vonsovskii 2J . The corresponding one- impurity 
model was applied later for the explanation of the 
Kondo effect the periodic version of this model was 
used for the description of the anomalous properties 
of strongly-correlated f-electron systems, e.g., heavy- 
fermion compounds fj]. The strong-coupling ("double- 
exchange" ) limit of the s-d model was applied for treating 
the interaction between localized moments and current 
carriers in a narrow d-band in ferromagnetic semicon- 
ductors (e.g., chalcogenide spinels) @ and manganates 
10. 

The strong interaction between the itinerant and lo- 
calized degrees of freedom in the above-mentioned class 
of systems can lead to incoherent structures in the spec- 
tral functions, in particular, the gap in electronic spec- 
trum can arise due to local correlations. Such a situa- 
tion occurs e.g. in the s-d exchange model which shows 
a metal-insulator transition with increasing the interac- 
tion of conduction electrons with localized spins Q. In 
more than three dimensions the local d — > oo picture @ 
is expected to give a qualitatively correct description of 
the physics near the metal-insulator transition, while the 
situation in three dimensions may be more involved. In 



particular, it was proposed that the appreciable short- 
range magnetic order in d = 3 may lead to the magnetic 
spin splitting above the Curie temperature Tq [ToL fill ] . 
Such a splitting, which can not be described within the 
local theory, was presumably observed in strong itinerant 
magnets Fe and Ni [l2l |. 

At the same time, in the peculiar d — 2 (or quasi-two- 
dimensional) case which is physically relevant for the lay- 
ered magnetic systems with small interlayer hopping (lay- 
ered manganites, cuprates, and other compounds with 
the perovskite structure), the strong short-range mag- 
netic order takes place. It is especially pronounced in the 
so-called "renormalized classical" temperature regime 
which arises above the magnetically ordered ground state 
and is characterized by an exponentially large correlation 
length £ cx exp(^4/T) {A is some constant) Strong 
magnetic correlations in this regime result in the struc- 
ture of the electronic spectrum, which is similar to that 
in the ordered phase [3, EH, EH- As a result, a quasi- 
splitting of the Fermi surfaces or a pseudogap can oc- 
cur at low enough temperatures due to ferromagnetic 
(FM) 17] or antiferromagnetic (AFM) /charge density 
wave (CDW) fluctuations|18j. These features of spec- 
tral functions have clearly different nature as compared 
to the abovediscussed Mott gap, since the finite (d = 2) 
dimension is crucial for their appearance; they are also 
expected to be relevant for quasi-2D systems with small 
interlayer hopping above the magnetic transition temper- 
ature. 

Recently, the pseudogap structures have been observed 
in layered manganite compound Lai 9Sri «Mn9Q7p^|. 
These structures are present both above and below Tc 
and possibly originate from the CDW fluctuations 
The FM fluctuations however may be responsible for the 
part of the pronounced shift (~250 meV) of the spec- 
tral weight maxima off the Fermi level above the Curie 
temperature. 

The separation of conduction and localized magnetic 
degrees of freedom can be also "dynamic" and arise in 
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purely itinerant electronic systems due to strong short- 
range magnetic order. In particular, it was pointed 
out that even purely itinerant electronic systems in the 
regime of strong magnetic correlations can be described 
by an effective classical s-d model [16j. 

To treat the effect of collective magnetic degrees of 
freedom on the electronic properties and describe prop- 
erties of 2D systems near magnetic quantum phase tran- 
sitions, another, spin-fermion model model, was intro- 
duced [2(| [21], [22j ■ Although this model was originally 
proposed as a phcnomenological model for systems with 
strong antiferromagnetic (AFM) fluctuations, the sys- 
tems with strong ferromagnetic (FM) fluctuations can 
be treated within this model as well [23|, [24| . 

The comparison of the spectral properties of sys- 
tems where magnetic fluctuations are induced by local- 
moments (described within the s-d model) or collective 
magnetic fluctuations (described within the spin-fermion 
model) in the regime of strong magnetic fluctuations is 
of certain interest. The spectral properties of the 2D 
systems with strong magnetic fluctuations were investi- 
gated previously within the spin-fermion model in the 
Eliashberg approximation [12, HII which neglects vertex 
corrections and in the quasistatic approach, performing 
summation of all the types of diagrams, but neglecting 
the effect of dynamic spin fluctuations [U [24[ ■ The lat- 
ter approach yields the two-peak structure of the spectral 
function for both strong FM and AFM correlations with 
a finite spectral weight at the Fermi level at finite £. 

Contrary to the spin-fermion model, where the length 
of the spin vector describing magnetic degrees of freedom 
at each lattice site is not fixed, the s-d model consid- 
ers fixed spin length. The result for the spectral weight 
A(kp,i/) ~ v 2 in the zero-bandwidth limit of the spin- 
fermion model [U [24[ shows that the spectral weight 
vanishes at v = only, and a true metal-insulator transi- 
tion is not expected to occur. Such a transition is however 
expected to be described by the s-d model [H, [2(| [13] • 
Treating of electronic properties near this transition is 
reminiscent the problem of considering the effect of mag- 
netic fluctuations in strongly-correlated Hubbard model. 
Although a number of approaches were proposed recently 
which extend the d — > 00 dynamical mean-field theory 
(DMFT) [Ito account for the short-range [H| and the 
long-range |29l. l30l . l3ll |32| correlations, the development 
of analytical approaches to this problem is also of certain 
interest. 

To study the magnetic fluctuations in the regime of 
the intermediate and strong interaction between localized 
and itinerant subsystems, the 1/z expansion, where z is 
the coordination number, can be used. This approach 
has been applied before to the ferro- and paramagnetic 

s-d model (see Refs. @, [M El E3) as wel1 as to the 
paramagnetic state in the Hubbard model [HI, (33[. Its 
extension to describe the paramagnetic state with strong 
magnetic correlations is of certain interest. 

Another promising candidate for non-perturbative de- 
scription of strongly-correlated electronic systems is the 



Ward- identity approach ^34]. Although this approach was 
quite successful in describing the strongly FM ordered 
state UK , its generalization to the weak FM and param- 
agnetic situation meets difficulties, since the contribution 
of longitudinal spin fluctuations, which should be taken 
into account in these cases, makes the system of equations 
for the electronic self-energy and electron- (para) magnon 
vertices not closed. 

This difficulty was overcome in a recently proposed 
combination of the Ward-identity approach with the 
1/M-expansion 17] (M is the number of spin compo- 
nents, M = 3 for the s-d and Hubbard model). The use 
of the 1/M expansion gives a possibility to truncate the 
hierarchy of integral equations for the self-energy and 
electron-magnon vertices. To apply this method to a 
broader class of models, we consider in the present pa- 
per its generalization - the equation of motion approach 
combined with the 1 /M expansion. Contrary to the 1/z- 
expansion, which considers equations of motion for Hub- 
bard operators, this method is applied to fcrmionic oper- 
ators, which allows us to describe both weak- and strong- 
coupling regimes, as well as the crossover between them. 

The plan of the paper is the following. In Sect. II we 
present the theoretical models describing interaction of 
electronic and magnetic degrees of freedom in strongly- 
correlated systems. In Sect. Ill we formulate the 1/z 
expansion for the s-d model. In Sect. IV we use the 
equation of motion method to derive the self-consistent 
equations for the self-energy and electron-magnon ver- 
tices within the 1 /M expansion, the details of the deriva- 
tion being presented in the Appendix. In Sect. V we 
investigate the electronic self-energy, spectral functions 
and electron-magnon vertices in different regimes. Fi- 
nally, in Sect. VI we summarize the main results of the 
paper. 

II. THE MODEL 

We consider a correlated electronic system with itiner- 
ant and localized-moment subsystems, which is described 
by the generating functional 

Z[rj\ = J D[c,J]D[S]exp{-S[c,S]/T 

-( c Lvka +vL c ka)} (1) 

where the fields c and S correspond to the electronic and 
spin degrees of freedom, respectively. The action <S[c, S] 
has the form 

l/T 

1 k 

q kk' aa' 

where first three terms describe the electronic and spin 
subsystems, the last term proportional to J corresponds 
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to the interaction between them; q = (q, ito n ), k = 
(k,ii/ n ); uj n — 2nnT and u n = (2n + l)irT are bosonic 
and fermionic Matsubara frequencies, £k is the electronic 
spectrum, <y aa i are the Pauli matrices, the term propor- 
tional to lZ q describes an indirect spin-spin exchange. 
Below we suppose that this exchange is ferromagnetic, 
i.e. 7Zq,o has its minimum at q = 0. 

Two different versions of the model © are considered 
in the paper. The first version, which is referred to as 
"s-d type models" has the measure of integration over 
the field S 



D[S] - £> sd [S] = ]JS(S 2 - S 2 )d 3 Si 



(3) 



so that the length of the field S is fixed to S. This is, in 
particular, the case for the lattice analog of the standard 
s-d exchange model 0] (i.e. Kondo-lattice model). In this 
case the vector potential A(S) describes the precession 
of free spin and satisfies the equation Vs x A(S) • S = 
1, Vs = (d/dS x ,d/dS y ,d/dS z ). In the classical limit 
S — ► oo this term fixes the field S to be static, dS/dr = 0, 
cf. Ref. (36|. For the s-d type of models the term pro- 
portional to lZq in Eq. ([2]) can be considered as arising, 
e.g., from the indirect spin-spin RKKY exchange. 

The model © with D[S] = D sd [S] and A(S) = 
(which is referred to as an effective s-d model in the 
following) can be naturally obtained from the Hubbard 
model in the renormalized classical regime [16]. The 
quantity lZ q is determined in this mapping by the RPA 
(random-phase approximation) expression for the inverse 
susceptibility in the ordered phase [l6| . 

Another version of the model © we consider is the 
spin-fermion model[20j, HJ, [22j with 



D[S] = Y[d 3 S u A(S) = 



(4) 



Contrary to the s-d model, the field S describes in this 
case three independent Bose-like fields, i.e. correspond- 
ing operators with different a commute with each 
other. The spin-fermion model can be also considered 
as that representing the low-energy electron and spin ex- 
citations of the Hubbard model in the quantum-critical 
regime [22l \3l\ . The quantity lZ q corresponds in this 
case to the inverse susceptibility Xq 1 an d enters as a 
phenomenological input parameter of the spin-fermion 
model. 

Both the static case dS/dr = and the case A(S) = 
allow the generalization of the model © to M- 
component spins S, = (Sj.-.S^ 1 ). This generalization 
will be used throughout the paper to perform the 1/M- 
expansion of the self-energy and electron-magnon ver- 
tices. 

For purely static lZ q — <5„o7?.q the s-d and spin-fermion 
models are described by the Hamiltonian 



H = ^ gkCko-Cka + ^ fcqSgS-q + H int) 



Hmt — —J ^ S 



k — k' G a a' ^ker ^k'er 



(5) 



Here the Fourier transformed operators Si which corre- 
spond to the fields in the continuum integral formal- 
ism, obey the standard SU(2) commutation relations and 
the condition S 2 = S(S + 1) for the s-d model, while the 
components of S commute with each other and S 2 is not 
fixed for the spin-fermion model. In the present paper we 
discuss only the classical limit of the s-d model with com- 
muting S-operators, so that the difference between these 
two models is only in the additional restriction for the 
S 2 for the s-d model (the quantum s-d model with the 
non-commuting spin operators will be considered else- 
where). For definiteness we assume J > 0; the results in 
the classical limit are the same for positive and negative 
J. 



III. THE 1/z EXPANSION FOR THE s-d MODEL 

The 1/z expansion considers perturbation theory 
around the atomic limit of the model (£k = const). 
The corresponding sequence of equations of motion for 
the electronic Green's function was considered in Refs. 
[2H [26j . To zeroth order in 1/z we obtain the result of 
the Hubbard-I approximation 



GkH 



1 



Fq{v) - £ k 



(6) 



where F (v) = v — I 2 /v, I = JS. To first order in 1/z 
one has to replace Fq(v) — * F-^(v) where 



F k (y) = v 



J 2 MT 

1 + E 



Xq 



1' ^ J F (^)-£k+ q 



(7) 



and x q is the non-uniform spin susceptibility of the model 
([5]). Despite the result (0 is obtained within the strong- 
coupling expansion of the s-d model, it reproduces cor- 
rectly in the second order in J the corresponding pertur- 
bation theory result for the electronic self-energy 0. 

The result ([7]) does not guarantee, however, correct 
analytical properties of the electronic Green's function at 
arbitrary J. To obtain analytical results, we consider self- 
consistent approximation of Hubbard-III type [25|, [26[ 

= v l + [F^)-F^)]G^) 



where 



1 + [y - F k (i/)]G k (i/) 



S 2 ^F k (^)-£ k+q ' 



(9) 



The solution of the Eq. © yields 



kk' 



l-^l + 4/ 2 G 2 H 
2 G k (v) 



(10) 
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The result for the electronic self-energy Sk(^) = v — 
F k (v), obtained from the Eq. flTDJ) with x q = S 2 /T, 
coincides with the result of the solution of single- impurity 
model of the dynamical mean-field theory (DMFT) for 
the classical s-d model, the Eq. ((9]) serving then as a 
self-consistency condition. 

IV. THE SELF-CONSISTENT EQUATIONS FOR 
THE SELF-ENERGY AND VERTICES WITHIN 
THE l/M EXPANSION 

To find the electronic self-energy in the model ([5|) 
within the l/M expansion, we consider the equation of 
motion for the fermionic operator c ka {r) in the Heisen- 
berg representation (see, e.g. Ref. [38| ) 

(d/dT + e w )d Wa (T) = [fltat.QcKr) (11) 

kV 

Applying this equation to the electronic Green's function, 
we find the expression for the electronic self-energy 

s^jEK + Vh, ( 12 ) 

where we have used the notation 

E'-^EE- 

q q uo n 

The vertices T describe scattering of an itinerant electron 
by the magnetic excitation (cf. Ref. [33|). 

T% q = J dT l dT 2 e ll '^ Tl ' T2)+lu - T2 (13) 

x(rpL(nFk +q .(r 2 )5 m q (0)])G,T 1 (G° fe+g )- 1 , 

T[...] stands for the imaginary-time ordering, G k and G k 
are the bare and full fermionic Green's functions, respec- 
tively, m = +, — , or z. Note that, despite the similarity 
of the definition of the vertex (fT3")) with that used for 
the Hubbard model in Ref. [34| . S*™ denotes in our case 
the spin-subsystem operator rather than the spin-density 
operator of itinerant electrons. 

The vertices (|13p can be found by applying equation 
of motion (jlip once more. For the following considera- 
tion, it is convenient to introduce the corresponding one- 
particle irreducible vertex with amputated bosonic and 
fermionic leg (cf. Ref. [341 ]) 

i£ n = Gl +q TlJ{J X ™G k+q ) (14) 

where x™ is the dynamic spin susceptibility. The equa- 
tion of motion relates the vertices (fl"4|) to the one-particle 
irreducible vertices of the interaction of the electron with 
two (para-) magnons. This vertex is in turn related to the 
vertex of an electron interaction with three paramagnons 



and so on (see Appendix for additional details). In the 
following we denote the vertex of an interaction of an elec- 
tron with n paramagnons as 7™^ "™" where mj = +, — , z 
are the paramagnon spin indices, k is the vector of the 
incoming fermion momentum and frequency, q\...q n are 
the vectors of the paramagnon momenta and frequencies. 
To truncate the resulting hierarchy of equations for the 
vertices, we follow the approach of Ref. [l7| and retain 
only the terms which contribute to the self-energy in the 
zeroth and first order in l/M (see Appendix). 

The resulting system of equations for the self-energy 
and vertices reads 

S fe = Mj2'lk +g G k+qXg (15) 
i 

7fe - l+Y.'i^-Mhl+qGl + q+lllqGk + ^Xq 
Q 

7 f = MY,'[H +q Gl +q + lk + qll z +q G 2 k+q 

+Jk Z +qG k + q ] Xq + UlkG k 

7fe 2± = -2M22' [lt +q Gl +q + ll +q ll z +q Gl +q 
i 

+lk+ql Z k + A q Gl +q ] Xq ~ UjlGl, 

where G k = [iu n — £k — £fc] _1 is the full (dressed) elec- 
tronic Green's function, 7™ 1 '" m " = 7™q "™" are the 
electron-paramagnon vertices at zero paramagnon mo- 
menta, 7fc = 7|. = 7^ , Xq = Xq, a = for the sprn- 
fermion model and a = — 2 for the s-d model. 

Comparing Eqs. (fl~5|l with the result (fT0|) of the 
1/z expansion we see that these expansions treat ver- 
tex corrections in different ways. While the l/M ex- 
pansion treats accurately the momentum dependence of 
the vertices, the "average" vertex of the 1/z expansion 

[\Jl + 4/ 2 Gj^)-l]/[2/ 2 G£H] interpolates between the 
weak- and strong coupling regimes. Both approaches, 
however, are expected to be applicable in the weak- and 
intermediate coupling regime. Below we consider the re- 
sult of application of these approaches to calculation of 
spectral functions in the s-d and spin-fermion models. 



V. THE RESULTS FOR THE SPECTRAL 
FUNCTIONS, SELF-ENERGY, AND VERTICES 

Equations (fl~5|) give a possibility to investigate the evo- 
lution of the spectral properties with varying electron- 
spin coupling J and the strength of magnetic correlations. 
For practical calculations, one has to specify an explicit 
form of the magnetic susceptibility Xq- I 11 the follow- 
ing we consider the low-temperature regime with strong 
magnetic fluctuations, where we employ an ansatz 

1 A 

x, = m?TP 4 "' |q|<<1 (16) 
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for the classical s-d exchange model and 



.4 



Xq = 



M q 2 +e- 2 +rK|/|q| 



lv F < |q| < 1 



for the effective s-d and the spin-fcrmion models (we have 
picked out explicitly the factor 1/M for further conve- 
nience). While the constant A is determined by the stiff- 
ness of spin excitations and can be arbitrary, the corre- 
lation length £ for the s-d type models is chosen to fulfill 
the sum rule 



MT^X q = S 2 



(18) 



For the spin-fermion model the correlation length is an 
independent parameter, which, however, can be related 
to the quantity A according to 



MT]r X(q , 0) = (A/j) 2 



(19) 



For the ferromagnetically ordered ground state A is ex- 
pected to be almost temperature independent at low 
T . We also consider the high-temperature regime with 
Xq = S 2 /(MT)S n0 . 

Note that the ansatz (fT6|) and (fT7|) neglects the non- 
analytical corrections to the spin susceptibility [4(| which 
are expected to be not too important in the renormalizcd 
classical regime. 



A. T = and high-temperature results 

We start investigating the solutions of Eqs. ([9]) and 
(TT5)) from the £ — * oo (T — > 0) limit. In this case the 
dominant contribution to momenta sums comes from the 
vicinity of q = point. Neglecting the momentum and 
frequency transfer q in the Green's functions and the ver- 
tices in the Eqs. ([9|) and ifPojl we obtain for the s-d model 
in the both 1/z and 1/M expansion, 



j2 j2 
Sfe = — , 7fc = 1 - =2 



(20) 



where V = v — + iO. The result (|20|) with V — v 
coincides with that in the atomic limit of the classical s-d 
model, where the spin dynamics can be indeed neglected 
and there are two energy levels at v = ±1 for parallel and 
antiparallel projections of the spins of itinerant electron 
and localized moment. In the £ — > oo case the two poles 
of the Green's function correspond to the splitting of the 
Fermi surface by strong magnetic fluctuations. 

The result ([2H)l can be compared to the corresponding 
result of the 1 /M expansion for the spin-fermion model 



St. — 



Ik 



M(Ag 



- ^v 2 - ai A 2 ^V 2 -a 2 A 2 ) 



2(2 + M) V 



M 



2—— [27T 4 + (6 + M)Ap 2 

2 (2 + M) A 2 V 2 1 



-MAq + (MA 2 , - 2V 2 



a x A 2 



(21) 



where A = A(T -► 0), a h2 = 1 + 4(1 ± y/l + M/2)/M; 
the branch Imy^z > of the square roots is chosen to 
guarantee the correct analytical properties for S and 7. 
The result (|2"Tj) yields also the two-peak structure of the 
spectral function with peaks at v ~ ±Ao and a finite 
gap A s f ~ A /2 at the Fermi level. As discussed in 
Refs. [It], the finite gap in the spectral function 
in this case is however an artifact of the first-order ap- 
proximation in 1/M, since the actual spectral function 
in the atomic limit has a behavior at small frequencies 
A(v) ~ |z/| M_1 which is non-analytic in 1/M (see Refs. 
|2ll[24|). However, as it will be shown below, the gap dis- 
appears quickly at finite £, making the situation in this 
case more favorable for the application of 1/M expansion. 

Now we study the solutions to the equations (TrS)) at 
finite £ (and T). At high temperatures the correlation 
length is small and the q-dependence of the magnetic 
susceptibility can be neglected, \q — S 2 /(MT)S n Q. The 
criterion for this behavior is T 3> J 2 /vf for the s-d ex- 
change model and T 3> vp for the spin-fermion model 
(vf is the Fermi velocity). The latter criterion is never 
fulfilled in real systems (which is a consequence of the 
fact that the local moments are not formed in the spin- 
fermion model, i.e. the magnetism formation mechanism 
is essentially Stoner-like), so that we consider the s-d 
model only. In the abovementioned limit we obtain in 
the 1/z expansion 



- [ v /l + 4/2(Gioc)2_ 1]/(2G! ioc ) (22a) 

/OO 
p Q {v')(v-v' (22b) 
-00 

where po{v) is the non- interacting density of states for the 
dispersion e^. As mentioned in Sect. Ill, Eq. (I22ap coin- 
cides with the result of the solution of the corresponding 
impurity model of DMFT, while Eq. (|22bj) is the self- 
consistency condition. The result (|22[) is correctly repro- 
duced by the 1/M expansion in the d — > 00 limit which 
is obtained by replacing the non-local Green's functions 
by the local one, G l ° c = X)k^ fe - At finite d this expan- 
sion leads, however, to different result, which we do not 
consider in this paper. 

The evolution of the self-energy £„, the spectral func- 
tion at the Fermi surface A(v) = — Im G^.p/ix, and the 
density of states p(y) — — Im G l ° c /it with increasing I 
calculated according to Eq. (|2"2"|) for the Bcthe (semiel- 
liptic) bare density of states poiy) — 1\J D 2 — v 2 /(irD 2 ) 





FIG. 1: The spectral functions (a), the interacting density 
of states p{y) (b), and real and imaginary parts of the self- 
energy (c,d) of the s-d model in the high-temperature regime 
at I = O.lnD, n = 1...5 (a,b) (from below to above), n = 2. ..7 
(c) and n = 1...6 (d) (from above to below). 



is shown in Fig. 1 (D is half of the bandwidth). Although 
the spectral functions have one-peak structure at small 
enough /, one can observe that the quasiparticle picture 
is in fact invalid at arbitrarily small / since the real part 
of the self-energy has positive slope and |ImE„| is max- 
imum at the Fermi level. These features are very sim- 
ilar to those observed earlier in the spin-fermion model 
in the presence of strong magnetic correlations 24j, al- 
though their physical origin in the present case is dif- 
ferent. Above a critical value I > I c = 0.5D both the 
real and imaginary parts of the self-energy diverge at the 
Fermi level and the spectral functions have the two-peak 
structure. The density of states p{v) acquires a gap at 
the Fermi level above I c , which corresponds to a metal- 
insulator transition. 



B. Low-temperature results in the weak-coupling 
regime 

Now we investigate the low-temperature regime with 
large correlation length £ and the susceptibility ansatz 
Q160 or (|17p. First we consider the contribution of the 
static spin fluctuations with zero bosonic Matsubara fre- 
quencies ui n = only (so-called static approximation). 
The effect of the dynamic spin fluctuations with non-zero 
bosonic Matsubara frequencies, which are less important 
in the renormalized classical regime, is considered below. 
In the static approximation and for the nearly constant 
Fermi velocity vp in the vicininity of the Fermi momen- 
tum kp both the self-energy and vertices are actually 
functions of iv n — only: 



H k = Y,{iv n - £ k ); 7 £ 



m\..m Tl m\..m 



7 ™--™»(«,„- £k ) (23) 



The integration over momenta can be then simplified by 
introducing an auxiliary variable a — upq(k F /£;p) (cf. 



FIG. 2: Spectral functions (a) and the imaginary part of the 
self-energy (b) of the classical s-d model in the weak-coupling 
regime (I — 0.1vp,D — 5vp) at different low temperatures 
within the l/z-expansion 



Ref. [24]), 



AT 



D 



2^7 fc+q ) ^/c + qXq J ^2 + „2 £-2 



da 



-D 

[ 7 mi - m "(i/-a)f 
[v - a - E(i/ - a) + i0] r 



(24) 



where we have performed analytical continuation iv n — > 
v + iO and restricted the integartion over a by ±D to 
account for the effect of finite band width. 

To solve Eqs. ©, (fTOjl . and (fl~5j) numerically we 
parametrize the energy dependence of the self-energy 
E(i/) and vertices -y m i-- m » (j/) by a set of values at the 
points {v r }r=i suitably chosen on the real axis (we choose 
N e ~ 800 -j- 1000) and iterate these equations with some 
initial condition until the convergence is reached. In the 
following we put A = 1/vf for defmiteness. 

First we consider the results of the solution of these 
equations for the s-d model. The frequency dependence 
of the self-energy T, u and the spectral functions within 
the l/z expansion is shown in Fig. 2, the correspond- 
ing results for the self-energy, specral functions and the 
vertex 7^ within the 1/M expansion are shown in Fig. 3. 

One can observe that the form of the spectral func- 
tions in l/z and 1/M expansions is similar. In particu- 
lar, we find that the imaginary part of the self-energy has 
a Lorentz-like form with |ImE(0)| ~ T£, the real part of 
the self-energy acquires a large positive slope at the Fermi 
level, <9ReE/<9^ ~ T£ 2 , so that the quasiparticle picture 
is invalid. With decreasing inverse correlation length to 
the values i;" 1 ~ Ao/vf, the structure of the spectral 
functions changes from the one-peak to two-peak form. 
At £ -1 < the electronic spectrum and correspond- 
ing Fermi-surface are quasi-split due to strong magnetic 
fluctuations. 

These results are similar to the previously obtained 
results for the spin-fermion model within the quasistatic 
approach, Ref. pj, H Hi]. The results of the 1/M- 
expansion for the spin-fermion model are shown in Fig. 
4. At large £ we obtain A = [{AT J 2 / Att) ]n(D£/v F )] 1/2 ; 
this quantity approaches a constant value Aq (which is 
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FIG. 3: The spectral functions (a), the imaginary and real 
parts of the self-energy (b,c), and the vertex function 7 (d) of 
the s-d model in the weak-coupling regime in the first order 
1/M-expansion at M = 3, I = O.Ivf, D = 5v F and different 
values of teperature. 



the analog of the interaction I in the s-d model) in the 
T — > limit in the renormalized-classical regime. The 
qualitative behavior of the self-energy, spectral functions 
and vertices is very similar to those of Fig. 3 for the s-d 
model and the results of Ref. [24[. In particular, the 
spectral functions have two-peak form at small enough 
£ _1 with the peaks being, however, wider than found for 
the s-d model. 

Now we consider the effect of the dynamic spin fluctu- 
ations with nonzero bosonic Matsubara frequencies. In 
the approximation of constant Fermi velocity v F the self- 
energy and vertices depend on the difference iu n — £k- 
After calculation of the self-energy at the imaginary- 
frequency axis we use the Pade approximants to perform 
analytical continuation of the results to the real axis [3^ | . 
Alternatively, one can perform calculations directly on 
the real axis, the results of both approaches being al- 
most identical. The integration over momenta in Eqs. 
(TTSl) can be performed similar to Eq. ([24]) . 



k+qXq 



q ILJ„ 



AT 
4ttM 



d 



-D 



/ M „(a)[7 mi - m "("'n+^n-fl)] P 

iu n + iu) n — a — E(ii/ n + iuj n — a)]' 1 



D/vf 



qdq 



A 



* J vV - (a/v F ) 2 q 2 + £- 2 + r\u n \/q 

\a\/v F 



(25) 



The results for the self-energy and spectral functions 
in the renormalized classical regime of the effective s- 
d model with different bandwidth are shown in Fig. 
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FIG. 4: The same as Fig. 3 for the spin-fermion model with 
J = 0.2v F , A = 0.1v F . 



5. Although the low-energy behavior of the self-energy 
and spectral functions (Figs. 5a, b) is similar to that 
found in the static approximation (Fig. 3), the self- 
energy at higher frequencies Ao -C H D behaves as 
Ti{y) ~ v 2 l z . At very high frequencies \v\ > D the behav- 
ior of the self-energy S(i/) ~ 1/v is restored. The results 
for the spin-fermion model with the dynamic spin fluc- 
tuations are shown in Fig. 6 and have the same qualita- 
tive frequency dependence as for the effective s-d model; 
the low-energy behavior of the spectral functions (Figs. 
6c, d) is close to the corresponding results of the static 
approximation, Fig. 4. Therefore, the effect of dynamic 
spin fluctuations in the renormalized classical regime is 
not pronounced, and the static approximation describes 
correctly the low-energy behavior of the self-energy and 
spectral functions. 

A treatment of the dynamic spin fluctuations with 
nonzero bosonic Matsubara frequencies is necessary to 
consider the self-energy and spectral functions in the 
quantum critical regime. Although in that case the 
ansatz for the susceptibility (|16p may not be correct be- 
cause of the possible non-analytic corrections [40|, it is 
still interesting to investigate the self-energy and spec- 
tral functions in this regime with this ansatz to clarify the 
role of dynamic spin fluctuations. In the quantum-critical 
regime the quantity A = [{AT J 2 /4vr) \n(D£/v F )} 1/2 be- 
comes temperature-dependent itself and one has to em- 
ploy additional ansatz for the temperature dependence 
of correlation length. According to the Hertz-Millis the- 
ory [4l| (which is valid in the absence of non-analytical 
corrections to susceptibility), 



r 



B{T/v F ) x '\ 



(26) 



for definitness we put B = 1. 

The resulting self-energy and spectral functions are 
shown in Fig. 7. The frequency dependence of these 
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FIG. 5: The spectral functions (a), the imaginary and real 
parts of the self-energy (b-d) of the effective s-d model in the 
renormalized classical regime with account of the dynamic 
spin fluctuations at the bandwidth D — vf (a-c) and D = 5vf 
(d), / = O.Ivf, and different values of temperature. 
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FIG. 7: The spectral functions (a,c) and the imaginary and 
real parts of the self-energy (b,d) of the spin-fermion model 
in the quantum-critical regime for different values of the cor- 
relation length. 
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FIG. 6: The same as Fig. [5] for the spin-fermion model with 
J = 0.2v F , A = O.Ivf. 



quantities is determined by the interplay of static and 
dynamic fluctuations. At small the dynamic fluctua- 
tions dominate, the real part of the self-energy has nega- 
tive slope and the imaginary part has a minimum at the 
Fermi level, the spectral functions having one-peak struc- 
ture, so that the quasiparticle picture is valid. However, 
the behavior of the self-energy ~ v 2 ^ at interme- 
diate frequencies (Tvp) 1 ^ 2 <C \v\ <C D leads to strong 
damping of the quasiparticles, similar to the charge in- 



stability case [42]. At larger the static fluctuations 
start to dominate and the low-frequency behavior found 
in Ref. [24| in the quasistatic approximation is restored: 
the imaginary part of the self-energy acquires maximum 
at the Fermi level, the real part having a positive slope; 
this invalidates completely the quasiparticle picture. Al- 
though the two-peak structure of the spectral functions 
is formed at intermediate only a structure with one 
broad peak survives with increasing The observed 
behavior of the self-energy is also very similar to that 
found for the Pomeranchuk instability case in Ref. [i^ ]. 
This shows that the vertex corrections (which are differ- 
ent in the cases of charge and spin instabilities) are not 
too important in the quantum critical regime. 



C. Low-temperature results in the 
intermediate-coupling regime 

With increasing the interaction / in the s-d model one 
approaches the metal-insulator transition. To obtain nu- 
merical solutions in this regime, it appears necessary to 
use smooth cutoff corresponding to the finite bandwidth, 
e.g. introducing the cutoff function 



w[a) 



tanh 



2\a-D\ 
D 



(27) 



in the integration over a in the Eq. (|2"4"]) and extend- 
ing the limits of integration from — oo to oo. The fre- 
quency dependence of the spectral functions A{v) and 
the self-energy £„, obtained from the l/z expansion for 
D/vp = 2, 1,0.7 is shown in Fig. 8. One can see that 
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FIG. 8: The spectral functions of the classical s-d model in the 
intermediate-coupling regime I = 0.5vf for the bandwidth 
D — 2vf (a), D — tip (b), and D = Q.7vf (c) at different low 
temperatures, (d) shows the imaginary part of the self-energy 
for D — vf- 



at temperatures lower than some finite temperature T s 
(i.e. at large enough £) the spectral functions vanish at 
small is, which corresponds to the splitting of the elec- 
tronic spectrum by magnetic fluctuations. The finiteness 
of the temperature T s is, however, most likely an artefact 
of l/z expansion. The imaginary part of the self-energy 
also vanihsh at small v and T < T s , except for the point 
v = 0, where it has a 8- function contribution — a8(is), a 
being some constant. 

The results of the 1/M expansion for the spectral func- 
tions A(y) and the self-energy E„, and vertices 7„, at the 
Fermi surface at M = 3 are shown in Fig. 9. One can ob- 
serve that at not too low temperatures the results of 1 /M 
expansion are similar to the results of l/z expansion. At 
low temperatures the 1/M expansion shows finite spec- 
tral weight at the Fermi level, but has additional spikes of 
ImE(^) at \ v\ — v c — Ao, at some temperature T' s we find 
ImS^c) = 0, Re£(zy c ) = v c , so that the pole of the elec- 
tronic Green's function shifts to the real axis. At T < T' s 
the solution to the equations (fT3|) becomes non-analytical 
in the upper half-plane. Although small damping of elec- 
tronic excitations is expected to appear in higher orders 
of the 1/M-expansion, the temperature T' s seems to cor- 
respond to a crossover to the regime with strongly coher- 
ent excitations at new preformed Fermi surfaces. 



VI. CONCLUSION 

We have investigated the problem of splitting of the 
electronic spectrum due to coupling of electronic degrees 




FIG. 9: The spectral functions of the classical s-d model with 
the bandwidth D — 2vf (a) and D — vf (b) in the first order 
1/M-expansion at M — 3 at different low temperatures and 
/ = 0.5vf- (c-f) show the real and imaginary parts of the 
self-energy and the vertex function 7 for D = 2vf 



of freedom to local magnertic moments and collective 
magnetic modes in itinerant systems. In both regimes 
of weak- and intermediate coupling of magnetic and elec- 
tronic degrees of freedom the spectral functions at suf- 
ficiently low temperatures have a two-peak form, aris- 
ing as a pre-splitting of the energy spectrum. In the 
intermediate-coupling regime, the width of the peaks in 
the classical s-d model tends to zero at the temperature 
T' s in the first order of the 1 /M-expansion, making the 
excitations at the preformed Fermi surfaces coherent at 
T < T' s . Although a small damping of the electronic 
excitations is expected to arise at finite T in higher or- 
ders of the 1/M-expansion, it is not expected to change 
qualitatively the form of the spectral functions. In l/z 
expansion we observe the full splitting of electronic spec- 
tra in the intermediate coupling regime at sufficiently low 
temperatures. 

The dynamic spin fluctuations do not change qual- 
itatively the low-frequency behavior of the spectral 
functions in the renormalized-classical regime at weak- 
coupling. These fluctuations are nevertheless important 
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in the quantum-critical regime, where they lead to a rich 
behavior of the spectral functions, which is determined 
by the intetrplay of static and dynamic fluctuations. 

A generalization of the results obtained in the present 
paper to the quantum s-d model is of great interest. This 
generalization will allow to describe analytically the for- 
mation of the Kondo resonance at the Fermi level near 
the metal-insulator transition in d — > oo, as well as 
its interplay with magnetic degrees of freedom in finite- 
dimensional quantum Kondo lattices. Another problem, 
which is important for the theory of correlated magnetic 
metals, is the investigation of the local-moment forma- 
tion and magnetic fluctuations in itinerant-electron sys- 
tems (e.g., described by the Hubbard model). 

We are grateful to P. Igoshev for careful reading of 
the manuscript. The research described was supported 
in part by Grant No. 4640.2006.2 (Support of Scientific 
Schools) and 07-02-0 1264a from the Russian Basic Re- 
search Foundation. 
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FIG. 10: Diagram representation of Eqs. (J33J) . Solid lines 
correspond to the bare electronic Green functions, the black 
circles stand for the one-particle reducible vertices (|32[) . wig- 
gly lines denote external legs of these vertices. The dashed 



are the reducible vertices of the interaction of an electron 
with n magnons, 



Appendix. The derivation of the equations for the 
self-energy and vertices 



■"igi...g„ — 



J dT l ...dT n e iuiTl+ - +l ^ 



In this Appendix we consider the derivation of the sys- 
tem of integral equations for the self-energy and vertices. 
Using the equation of motion (jlip we obtain 



H (d x )(T[c x clSrSpST 3 -] = S(x - y) (28) 
x (T[SrsrST 3 -}) + J(T[c x 4(aS x )S^Srsr 3 ...}) 

where Ho(d x ) = d T + e(V), x,y,a,b, and c are the space- 
time coordinates. Translating this to the momentum- 
frequency space and combining the result with the Dyson 
equation G _1 = —(Ho + £) we find 



(29) 



pmi . . .m 71 ;erer 
k;qi...q n 



— RTOi...ro„ _ rp j 

^fc °qi+... + q n rL qi...q„ 1J 

■Hm\—m n+ i,aa" m n + 1 
1 fc;«i. ..<?„+! a a"a' 

q rl + 1 ,m Tl + 1 ,a" 
ALjr fc+9l + ...+9„+l 



where G° k = (iv n — £k) 1 , 



x(T[^(n)...^(r n )]) (31) 

is the spin correlation function. 

To obtain the equations for 1-particle irreducible ver- 
tices we first define the connected vertices via the recur- 
sive relations 



k;qi...q n ~ k;qi...q n 

■tOn-I 



(32) 



+(G k ) 5 qi+ ,.. +qn R'^ 1 ;; q " 



EEv 

s=2 {/»} 



h9i 3 ■""«»! — q* a 



xT k^ s + 1 ...q Zn 

where 7™ = {£ r }£_j is the s-element ordered subset of 
{l...n}; {ir}r=s+\ = {l--- n }\^s l j an< i t ne summation is 
taken over all subsets 7™. The equations for the vertices 
(1321) read 



JJ L \ R qi...q n .,qi+-+q n (J <?<?' 

m„ + i 

n-1 

—rriij ...m is ,m„ + i 

2 / n qi 1 ---qi s ,qi 1 +---+qi s 

s=l {I s },(t" 



pmi .. .m^jcrcr 
k;qi...q n 



K^TqT" = -G^(Gl +qi+ ... + J-i j dr.dn...dr n 



x(T[^(rx)...^(T n )4 ff (r) (30) 
q„,<T'(0)])e il/T+Miri+ -+^" r " 



xc k 



-«Ji- 



XG k+q zl +...+q Zs T k,q is + 1 ...q in 
L k;q 

971 + 1:0-" 

X G 1 + qi +. ..+<?„+! } 



/a i . +1 ...Tni„,o-CT Q.rnm+1 



<7i...<j„ + i a a"a' 



(33) 
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where R denotes the connected part of R. The diagram 
representation of the Eqs. for n = 1...3 and the Eq. 
([29| for the self-energy is shown in Fig. [lOl 

At the next step we express the vertices r™^"™" 
through the corresponding one-particle irreducible ver- 
tices 7™^"™" with the use of the Legendre transforma- 
tion of the generating functional, see, e.g., Ref. [43j]. To 
the first order in 1/M the equations for ^k, qi , lk Z qiq2 > 
and 7fc Z g f g2(J 3 form closed system, see Fig. 11. The ana- 
lytical form of these equations reads (below we omit the 
electronic spin indices of 7 assuming that any of the equal 
nonzero spin components of the vertex is taken): 



Sfc — M 'lk;qGk+q 



Xq 



(34) 




7fc,<?l = 1 + ' t( 2 ~ M hk;q 2 lk+q 2iqi G k+q2 

92 



9192 



~^~1k;qi ,92] X92^fc+gi+92 

M ^2 'Gk+q 1 +q 2 +q 3 Xq 3 



(35) 



FIG. 11: Diagram representation of the Eqs. (J34J) - <j3Sj) - White 
circles stand for the one-particle irreducible vertices 7™,j"™" , 
dashed circles for the vertex ^~^-[~ q2 q 3 , Eq. (|38[l . bold solid lines 
correspond to the full electronic Green function, wiggly lines 
- to the magnetic susceptibility \q- 



lk;q 1 q 2 q a 



X [7fc;<?3 lk+q 3 \qi lk+q 3 + qi ;q 2 Gk+q 3 Gk+q 3 + qi 

+7fc;g 3 7fe+93;92 7fe+93+«2;9i Gk+q 3 Gk+q 3 +q 2 

l ZZ± _r_ ZZ (~i 

" 1 ~lk,q 1 q 2 q 3 > lk,q 3 !k+q 3 ,gig 2 ^fc+gs 
+Mjk,q 1 +q 2 +q 3 Tq iq *q 3 Xq 1 +q 2 +q 3 ] (36) 

= Gk+q!+q 2 +q 3 +qiXq4 



<ll 



Ylk\qi {jk+q 4 ;q 1 1k+q 4 +q 1 ;q 3 Jk+q 4 +q 1 +q 3 ;q 2 

X Gk+q 4 Gk+q 4 + qi Gk+q 4 + qi +q 3 

"^"/k+q 4 ;q 2 r ^k+q 4 +q 2 ;q 3 "/k+q 4 +q 2 +q 3 ;q 1 

X Gk+q 4 Gk+q 4 +q 2 Gk+q 4 +q 2 +q 3 

~ lk+q 4 ;g 3 7fc+94 +93 191 7fc+94+gi +93 ;<?2 

X Gk+q 4 Gk+q 4 +q 3 Gk+q i + q3 +q 1 

— 7fc+g4;93 7fc+94 +93 ;92 7fc+94+92 +93:91 
x Gk+q 4 Gk+q 4 +q 3 Gk+q i +q 3 +q 1 

~lk+q 4 ;q 4 7^+94 +9l ;92 7fc+94+9l +92 ;93 
'XGk+q 4 Gk+q 4 + q 1 Gk 

+94+91+92 
— 7fe+94;92 7fc+94+92;9l7fc+g4+9l+92;93 
X Gk+q 4 Gk+q 4 +q 2 Gk+q 4 + qi +q 2 

~7fc+g 4 +9l +92 ;93 lk+q 4 ;q 1 q 2 Gk+q 4 Gk+q 4 + qi +q 2 
~lk+q 4 ;q 3 r yk+q 4 +q 3 ;q 1 q 2 Gk+q 4 Gk+q 4 +q 3 

~"fk+q 4 ;q 1 q 2 q 3 Gk+q 4 ) ~ M~/k+q 1 +q 2 +q 4 ;q 3 

X ^91 92 94 Gk+qi +q 2 +q 4 Xqi +92 +94 7fc,9l +92 +94 

— 7fe+9l +92 +94 ;93 7fe; 9l q 2 q 4 Gk+q 4 +g 2 +94 

(37) 

where Xg = Xq, 

j^mim 2 m 3 m 4 ti t— 2 rymirn 2 m 3 m 4 

gi9293 gi9293,gi+92+g3 

X (\ mi y m! v m3 Y m4 f38l 
X Ugi Ag 2 Xq 3 X Ql + q2 + q3 ) [J ) 

is the irreducible vertex of the two-magnon interaction. 
In previous treatments of the spin-fermion model (see, 



e.g., Ref. [2l|) the vertex ([38j) was neglected. For the 
flat density of states it indeed vanishes in the limit of 
zero bosonic frequencies. Although it was pointed out 
that this vertex is singular at nonzero external frequen- 
cies and T = this result is of possible relevance at 
finite temperatures in the quantum-critical regime only, 
since in the renormalized classical and high-temperature 
regimes the most important contribution to the self- 
energy and vertices comes from the terms with zero 
bosonic Matsubara frequencies. In this regimes the ver- 
tex (|38[) may be non-negligible for the non-constant (es- 
pecially singular) density of states. At the same time, 
the vertex ([38|) does not vanish even for the flat density 
of states in the s-d model. To find this vertex we apply 
the 1/M-expansion [H| to the spin part of the model ([5]). 
To zeroth order in 1 /M we find 



Y+ — ZZ _ n 

qiQ2q3 



(JMfTj^Xq'Xq'+q^ 



(39) 



To perform numerical calculations we consider two 
cases (see also the main text): (i) the susceptibilities 
X q are static and momentum-independent, Xq = &noXi 
the electronic Green's functions also do not depend on 
the momentum, Gk = Gw n (high-temperature regime); 
(ii) the susceptibilities have sharp maximum of the width 
at q = 0, the momentum and frequency dependence 
of the electronic Green's functions being arbitrary (low- 
temperature regime). 

In both the cases we approximate in the equation for 
the vertex 7™i"'™" (n = for the self-energy) the Green's 
functions and vertices at qi ^ by the corresponding 
quantities at the maximum total 4-momentum k + q\ + 
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... + q n , i.e. put 

mi...m r ^ mi...m r 

7 ^+Ei mQi-Qr ~~ '/c+ 9 i+. ..+g„ + i 

G fe+E 4 <?» — .+«„+! (40) 

In addition we use the approximation 

i^XQ+gsXift 1 q 1 q 2 q3 

93 93 

= -2/ 9ii o (41) 

where /g i; q 3 is some function of electronic Green's func- 
tions and vertices, and we have taken into account that 



^<jig2<23 depends in fact only on two first momenta. 

Eqs. (|40p and (f4"Tj) are exact in the regime (i) where the 
electronic Green's functions and vertices are momentum- 
independent and the bosonic Matsubara frequencies can 
be put to zero. In the regime (ii) the approximation (|40j) 
underestimates slightly the effect of magnetic correla- 
tions, but treats them qualitatively correct and becomes 
exact in the limit of infinite correlation length (where 
only vertices with the momentum q — enter). Approxi- 
mation (|41[) is even more accurate in the regime (ii), since 
the product XQ+qsXqs is even stronger peaked at q$ = 0. 
Applying (gDJl and (@TJ to ([51 )1 - ([57 ]) we obtain the system 
of equations (fH>)) of the main text. 
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